Observation of elastic anomalies driven by coexisting dynamical spin Jahn- Teller effect 
and dynamical molecular spin state in paramagnetic phase of the frustrated MgCr 2 04 
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Ultrasound velocity measurements of magnesium chromite spinel MgCr2 04 reveal elastic anoma- 
lies in the paramagnetic phase that are characterized as due to geometrical frustration. The temper- 
ature dependence of the tetragonal shear modulus (Cn — Ci2)/2 exhibits huge Curie-type softening, 
which should be the precursor to spin Jahn- Teller distortion in the antiferromagnetic phase. The 
trigonal shear modulus Caa exhibits nonmonotonic temperature dependence with a characteristic 
minimum at ^50 K, indicating a coupling of the lattice to dynamical molecular spin state. These re- 
sults strongly suggest the coexistence of dynamical spin Jahn- Teller effect and dynamical molecular 
spin state in the paramagnetic phase, which is compatible with the coexistence of magnetostructural 
order and dynamical molecular spin state in the antiferromagnetic phase. 

PACS numbers: 72.55,+s, 75.30.-m, 75.40.Gb, 75.50.Xx 
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I. INTRODUCTION 

Geometrically frustrated magnets have attracted con- 
siderable interest because the vast degeneracy of their 
ground states gives rise to a rich variety of novel 
phenomenal. One of the most highly frustrated spin sys- 
tems is the nearest-neighbor Heisenberg antiferromagnet 
on the pyrochlore lattice, where the magnetic sites form 
a network of corner-sharing tetrahedra. The prototypical 
examples for this highly frustrated system are chromite 
spinels ACr 2 04 with non-magnetic ions A = Mg, Zn, and 
Cd. The magnetic properties of these compounds are 
fully dominated by the Jahn- Teller (JT)-inactive Cr 3+ 
with spin 5 = 3/2 residing on the pyrochlore network. 
The high values of frustration parameter / = |0w/T/v|, 
the ratio of Weiss and Neel temperatures, indicate the 
presence of strong geometrical frustration: / ~ 30 with 
G w = -370 K and T N = 12.5 K for MgCr 2 4 , / ~ 33 
with G w = -390 K and T N = 12.0 K for ZnCr 2 4 , and 
/ ~ 9 with G w = -70 K and T N = 7.8 K for CdCr 2 4 2 . 

In ACr 2 04, the antiferromagnetic (AF) transition 
at T/v coincides with a tetragonal lattice distortion 
with a compression along the c-axis for MgCr 2 04 and 
ZnC^O^ 3 -^, and an elongation for CdCr 2 04— . This 
magnetostructural ordering is explained by the spin JT 
mechanism via the magnetoelastic coupling, where local 
distortions of the tetrahedra release the frustration in 
the nearest-neighbor AF interactions^. The very small 
magnitudes of the tetragonal lattice strain of the order 
10 -3 imply the presence of a nonuniform component in 
the lattice distortion, where the neighboring tetrahedra 
are alternately stretched and compressed along the same 
crystal axis 3 - — Indeed, a recent study of ZnCr 2 04 using 
synchrotron X-ray and neutron scattering experiments 
revealed complex spin-lattice order with the nonuniform 
lattice distortion^. 

In the paramagnetic (PM) phase of ZnCr 2 04 and 



MgCr 2 04, inelastic neutron scattering (INS) experi- 
ments observed quasielastic magnetic scattering, indi- 
cating the presence of strong spin fluctuations due to 
frustration 3 -^ - — . Interestingly, this quasigapless mode 
was characterized as fluctuations of AF hexagonal spin 
clusters (AF hexamers) in the pyrochlore lattic e) 11 ! 12 . It 
is natural to expect that this dynamical molecular spin 
state vanishes in the AF phase as a consequence of the 
release of frustration. Actually, however, the INS studies 
revealed the appearance of gapped molecular spin excita- 
tions in the AF phase, which indicates the coexistence of 
the magnetostructural order and the zero point motion- 
like spin moleculesii. 



The appearance of the magnetostructural order in 
ACr 2 04 means that the spin-lattice coupling plays a sig- 
nificant role in releasing the frustration; thus, the phonon 
spectra are expected to be affected by the frustration. 
Indeed, infrared and Raman spectroscopies in ZnCr 2 04 
and CdCr 2 04 revealed that the optical phonon spec- 
tra exhibit anomalies of softening on cooling in the PM 
phase, and splitting in the AF phasei 3 -— . In this paper, 
we present acoustic phonon information provided by ul- 
trasound velocity measurements in magnesium chromite 
spinel MgCr 2 04. The modified sound dispersions by 
magnetoelastic coupling allow one to extract detailed in- 
formation about the interplay between spin and lattice 
degrees of freedom in a lower energy region near to the 
ground states. Thus ultrasound velocity measurements 
can be a useful tool in studying geometrically frustrated 
magnets^"—. We find elastic-mode-dependent sound- 
velocity anomalies in MgCr 2 04 in the PM phase, similar 



to those observed in the past in ZnCr 2 04— . The ori- 
gins of these elastic anomalies are discussed in terms of 
frustration inherent in ACT2O4. 
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II. EXPERIMENTAL 

Ultrasound velocities were measured in single crystals 
of MgCr204 prepared by the floating zone method, where 
the phase comparison technique was used with longitu- 
dinal and transverse sound waves at a frequency of 30 
MHz. The ultrasounds were generated and detected by 
LiNbC>3 transducers glued on the parallel mirror sur- 
faces of the crystal. We measured sound velocities in all 
the symmetrically independent elastic moduli in the cu- 
bic crystal: compression modulus Cn, tetragonal shear 
modulus Cll ~ Cl2 = Ct, and trigonal shear modulus C44. 
The respective measurements of Cn, Ct, and C44 were 
performed by using longitudinal wave with propagation 
k||[100] and polarization u||[100], transverse wave with 
k||[110] and u||[110], and transverse wave with k||[110] 
and u||[001]. 
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III. RESULTS 

Figures l(a)-(c) depict the longitudinal sound velocity 
vl in Cn and the transverse sound velocity vt in C t and 
C44 as functions of temperature (T), respectively. All the 
elastic moduli exhibit a jump at T/v = 13 K. And, more 
importantly, all the elastic moduli exhibit elastic-mode- 
dependent anomalous T-dependence in the PM phase, 
from room temperature (300 K) down to T/v, that is dif- 
ferent from the T-dependence usually observed in solids, 
namely monotonic hardening with decreasing T— . These 
elastic anomalies should have magnetic origins where the 
spin degrees of freedom play significant roles but also 
taking into account the absence of orbital degeneracy in 
the 5-site Cr 3+ in MgC^C^. Such elastic anomalies are 
attributed to magnetoelastic coupling acting on the ex- 
change interactions. In this mechanism, the exchange 
striction arises from a modulation of the exchange inter- 
actions by ultrasound as follows >2i 



FIG. 1: (Color online). Temperature dependence of the lon- 
gitudinal (vl) and transverse (vt) sound velocities, (a) Cn, 

(b) (Cll ~ Cl2) = C t , and (c) C44. 



polarization u and propagation k relative to the exchange 
path 5. Concerning the elastic anomalies in MgC^C^, we 
emphasize that the anomalous T variations of the elas- 
tic moduli in the PM phase are divided into two types: 
Curie- type -1/T softening for Cn(T) and C t (T) and non- 
monotonic T-dependence with a characteristic minimum 
at —50 K for C^(T). Together these types imply that the 
elastic anomalies have different origins, which we shall 
discuss below. 



IV. DISCUSSION 
A. Curie-type softening in PM phase 



H exs = Y,[J(S + u l -u j )- J(5)}Si ■ Sj. (1) 

Here 5 = Ri — Rj is the distance between two magnetic 
ions, and U{ is the displacement vector for the ion Ri. 
When a sound wave with polarization u and propagation 
h is given by u = uoex.p[i(k • r — (Jt)], where uq and uo 
are the respective amplitude and frequency, the exchange 
striction of Eq. (pQ) is rewritten as>22 

f) j 

H exs = J2(qs ■ u o)(k ■ S)(Si ■ S^y^-^K (2) 

i 

Here the exponential is expanded to first order because 
with a 30-MHz ultrasound frequency kS <C 1. Eq. (j2j) 
implies that both the longitudinal and transverse sound 
waves can couple to the spin system via the exchange 
striction mechanism, which depends on the directions of 



First we discuss the origin of the Curie-type softening 
in Cn(T) and C t (T) in the PM phase. Taking into ac- 
count that the relative change in softening is much larger 
in C t (T) (^ -50%) than in C n (T) (^f -10%), 
and that the compression modulus C\\ is written as 
C\\ = Cb + f Ct with Cb the bulk modulus, the Curie- 
type softening should originate from a softening in the 
tetragonal shear modulus C t (T). Thus the Curie- type 
softening should be a precursor to the tetragonal struc- 
tural transition at T/v. Such a precursor softening to 
the structural transition is well known to occur as a re- 
sult of the Jahn- Teller (JT) effect in orbital-degenerate 
systems^. However, the absence of orbital degeneracy 
for the B-site Cr 3+ in MgC^CU rules out such orbital 
effects. Thus the Curie-type softening is most probably 
due to a spin JT effect 6,7 , which is compatible with the 
coincidence of AF ordering and tetragonal lattice com- 
pression at Tat. The spin JT distortion is analogous to 
orbital distortion in the sense that a degenerate electronic 
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system lowers its energy by lifting the degeneracy under 
structural distortion. In AC^O^ the spin JT distor- 
tion lifts the spin degeneracy (frustration) inherent in 
the pyrochlore lattice by modulating the balance of ex- 
change interactions via the exchange striction mechanism 
expressed as Eq. (pp). 

For the orbital JT effect, a JT- active ion and its set 
of ligands are considered to be a structural unit, where 
the crystal-field striction mechanism leads to the orbital 
JT distortion. T dependence of elastic modulus Cp(T) 
in the orbital JT systems is explained assuming the cou- 
pling of ultrasound to JT-active ions via the crystal-field 
striction mechanism, and the presence of inter- JT-active- 
ion interactions. A mean- field expression of Cr(T) in the 
orbital JT systems is given as^"— 



C r (T) 



9rN 



Xr(T) 



(l-5rXr(T))' 



(3) 



with Cp the elastic constant without the orbital JT ef- 
fect, gr the coupling constant of crystal field strain, N 
the number of JT-active ions in the unit volume, g r the 
inter- JT-active- ion interaction, and xr (T) the strain sus- 
ceptibility of a single JT-active ion. Eq. (j3j) expresses 
T dependence of the strain susceptibility which is the 
response function to an applied strain, and its form is 
analogous to the mean-field expression of the magnetic 
susceptibility which is the response function of the mag- 
netization to an applied magnetic field. In the orbital 
JT systems, the degenerate ground state is considered 
to couple strongly and selectively to the elastic modu- 
lus CV which has the same symmetry as the orbital JT 
distortion. For such a JT-active elastic mode, the strain 
susceptibility is dominated by Curie term Xr(T) ~ 1/T 
at low temperatures. Then Eq. (j3]) is rewritten as>22r— 



Cr(T) = C? 



T-T c 
T-0 



(4) 



with the inter- JT-active-ion interaction, and T c the 
second-order structural transition temperature. Here 
is positive/negative when the interaction is ferrodis- 
tortive/antiferrodistortive. Although the background Cp 1 
in Eq. (jl]) generally exhibits a hardening with decreasing 
T— , we here assume Cp 1 is constant because its hard- 
ening is negligibly small compared with the huge Curie- 
type softening in Cr(T). In contrast to the JT-active 
soft mode Cp expressed as Eq. (j4]), the elastic modulus 
Cr' with the different symmetry from the orbital JT dis- 
tortion is not affected by the degenerate ground state; 
Cp'(T) exhibits the absence of anomaly, Cr f (T) = Cp 1 ,, 
originating from xv = in Eq. ((3]). 

In analogy to the orbital JT effect, the Curie- type soft- 
ening in Ct(T) in MgC^CU can be explained in the spin- 
JT-effect-scenario assuming the coupling of a Cr 3+ struc- 
tural unit to the external strain of ultrasound via the ex- 
change striction mechanism expressed as Eqs. (pQ) and (j2j), 
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FIG. 2: (Color online). The tetragonal shear modulus Ct in 
MgCr2 04 (open circles, from Fig. 1(b)) and ZnC^CU (open 
squares, from Ref. [ 16 ]) as functions of temperature. The 
solid and dotted curves are fits of the experimental data for 
MgCr2 04 and ZnC^CU to Eq. (4]), respectively. Values for 
the fitting parameters are listed in the table beneath the fig- 
ure. The inset shows a schematic of the tetragonal lattice 
deformation in Ct- 



and the presence of inter-structural- unit interactions. Al- 
though the structural unit of the spin JT effect has not 
been identified so far, the mean- field expression of the 
soft mode Cr(T) should have the same form as Eq. 
with Cp 1 the elastic constant without the spin JT effect, 
the inter-structural-unit interaction, and T c the second- 
order structural transition temperature. For MgC^CU, it 
is expected that ultrasound in Ct generates the tetragonal 
exchange striction which promotes the tetragonal spin- 
lattice correlation, while the trigonal exchange striction 
is generated in C44. Thus, since the spin JT distortion in 
MgCr204 has the tetragonal symmetry 4 , the tetragonal 
shear modulus Ct(T) should exhibit Curie- type soften- 
ing as expressed in Eq. (j4]), whereas its absence in the 
trigonal shear modulus is expressed by C±±(T) = C44 in 
analogy to the orbital JT effect. 

It is noted that the Curie-type softening in the 
tetragonal shear modulus C t (T) is observed in not 
only MgCr204 but also another highly frustrated spinel 
ZnC^O^. Figure 2 compares the tetragonal shear mod- 
ulus Ct in MgCr204 (open circles, from Fig. 1(b)) with 
that in ZnC^CU (open squares, from Ref. f 16 ]) as func- 
tions of T; fits of the experimental data for MgC^CU 
and ZnCr204 to Eq. are also plotted as solid and 
dotted curves, respectively. The plots reproduce the ex- 
perimental data given the fitted values of C £ °, T c , and 
listed in the table beneath the figure. For MgC^CU and 
ZnCr204, the respective fitted values of T c = 0.7 K and 
K are lower than the experimentally-observed Neel tem- 
peratures of T/v = 13 K and 12.5 K, T c < T/v, indicating 
that the phase transition at T/v is of first order-* 4 *. The 
negative fitted values of mean the dominance of an- 
tiferrodistortive inter-structural-unit interactions which 
are stronger in MgC^CU (0 = -14 K) than in ZnC^CU 



(6 = -6 K). The dominance of such antiferrotype interac- 
tions in the PM phase seems to be compatible with the 
nonuniform lattice distortion in the AF phase revealed 
experimentally in ZnC^CU, where neighboring tetrahe- 
dra are alternately stretched and compressed along the 
same crystal axis^. In the antiferrotype nonuniform lat- 
tice order, the tetragonal lattice strain per spinel unit cell 
should be negligibly small, as is consistent with the very 
small magnitudes of the tetragonal lattice strain observed 
in ACr 2 4 in the AF phased-. 



B. Softening with minimum in PM phase 

Next, we discuss the origin of the other salient elastic 
anomaly in the PM phase, the characteristic minimum 
at ~50 K in C^(T). As shown in Fig. 3, this anomaly 
is observed in not only MgC^CU (open circles, from Fig. 
1(c)) but also ZnCr204 (open squares, from Ref. 
Taking into account the absence of orbital degeneracy 
for the 5-site Cr 3+ in these compounds, the spin degrees 
of freedom should play a significant role for the pres- 
ence of the anomaly in C^{T). One possible origin for 
the anomaly in C^(T) is the spin JT effect. However, 
the symmetry of the spin-lattice order in MgC^CU and 
ZnCr204 is tetragonal, which is different from the trig- 
onal symmetry of 644—^. Thus, in the spin-JT-effect 
scenario, as described above, while the tetragonal shear 
modulus C t (T) should exhibit Curie- type softening in the 
cubic phase (PM phase) as expressed in Eq. (J2J), the trig- 
onal shear modulus C44(T) should be anomaly- free as 
C44(T) = C44, which rules out the spin JT effect as a 
possible origin for the anomaly in C^(T). The other 
possible origin for the anomaly in C^(T) is the cou- 
pling between the dynamical molecular spin state and 
the acoustic phonons. We here note that the INS studies 
in MgCr204 as well as ZnC^CU observed the quasielastic 
magnetic scattering in the PM phase^i^— , and charac- 
terized this quasigapless mode as the fluctuations of the 
AF hexagonal spin clusters (the AF hexamers) in the py- 
rochlore lattic o 11 ! 12 . As illustrated in the inset to Fig. 3, 
the AF hexamers are aligned in the trigonal (111) planes 
of the spinel structure, and thus should be sensitive to 
the trigonal lattice deformation which generates the trig- 
onal exchange striction. Therefore, the anomaly of the 
trigonal shear modulus C±±(T) in the PM phase should 
stem from the dynamical molecular spin state. 

The ultrasound velocity measurements in the present 
study extract information about the acoustic phonon dis- 
persions at very small wave vectors of q ~ 10 -5 A -1 . In 
contrast, the molecular spin fluctuations observed in the 
INS spectra have large wave vectors of q ~ 1.5 A -1 — 
Thus the observation of the anomaly in C44(T) suggests 
that the strong coupling of the dynamical molecular spin 
state to the lattice modifies the acoustic phonon disper- 
sions over a wide range in the Brillouin zone, even at 
around zero wave vector. The modification of the phonon 
dispersions would be strongest at around q ~ 1.5 A -1 . 
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FIG. 3: (Color online). The trigonal shear modulus Caa in 
MgCr2 04 (open circles, from Fig. 1(c)) and ZnC^CU (open 
squares, from Ref. [ 16 ]) as functions of temperature. The 
solid and dotted curves are fits of the experimental data for 
MgCr2 04 and ZnCh^CU to Eq. ([5]), respectively. Values for 
the fitting parameters are listed in the table beneath the fig- 
ure. The inset depicts an AF hexamer in the Cr 3+ pyrochlore 
lattice. 



Regarding the elastic anomaly from the short-range 
spin correlations, the characteristic minimum in C^{T) 
in MgCr204 and ZnC^CU is similar to that in low- 
dimensional spin-dimer systems such as SrCu2(BOs)2r^. 
Notably, the minimum in Cr(T) in the spin dimer sys- 
tems is explained as a result of the presence of a finite 
gap for the local magnetic excitations which is sensitive 
to the strain^. In the spin-dimer systems, the quantita- 
tive analysis of the experimental data of Cr(T) is given 
by applying a theoretical model which has the same form 
as Eq. (j3j). Assuming a singlet-triplet excitation gap Ai 
for a single spin dimer and a multitriplet excitation gap 
A2, the contribution of the spin dimers to the elastic 
constant takes the formal 



C r (T) 



G{ 



2 N 



Xr(T) 



{l-K rXr {T)Y 



(5) 



with Cp the background elastic constant, N the density 
of spin dimers, G?i,r = |9Ai/<9er| the coupling constant 
for a single spin dimer measuring the strain (er) depen- 
dence of the excitation gap Ai, K-p the (q = 0) dimer- 
dimer interaction, and Xr(T) the strain susceptibility of 
a single spin dimer. Eq. (j5j) is based on the coupling 
between the gapped spin triplet states and the acoustic 
phonons, which, in turn, means that this equation can be 
applied to all systems with the coupling between gapped 
states and acoustic phonons regardless of whether it is 
a singlet system or not. According to Eq. (j5j), the mini- 
mum in Cr(T) appears when this elastic mode strongly 
couples to the excited state at A^; on cooling, Cr(T) ex- 
hibits softening roughly down to T ~ A^ , but recovery of 
the elasticity (hardening) roughly below T ~ A^. 

In analogy to the spin dimer systems in the presence 
of local excitations, therefore, the minimum in C^{T) in 
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MgCr 2 04 and ZnCr 2 04 should arise from a gap in the 
molecular spin excitations which is sensitive to the strain. 
This interpretation helps to understand the INS results 
for MgCr204 and ZnCr 2 04. The broad quasielastic mag- 
netic scattering spectrum in the PM phase comprises 
gapless spin fluctuations which probably form embryo of 
the tetragonal spin-lattice order below T/v, and gapped 
spin excitations which are considerably smeared^^— . 
And the observation of distinct gapped excitations in the 
AF phase is attributed to the suppression of spin fluctu- 
ations by spin JT distortion ^ 1Q i n . 

Similar to the spin dimer systems^, we now give a 
quantitative analysis of the experimental data of C±±(T) 
in MgCr204 and ZnCr 2 04 assuming an excitation gap 
Ai for a single AF hexamer and a multi-AF-hexamer ex- 
citation gap A2. The contribution of the AF hexamers to 
the elastic constant should take the same form as Eq. (j5j) 
with 644 the background elastic constant, N the den- 
sity of AF hexamers, Gi,r 5 = |dAi/<9er 5 | the coupling 
constant for a single AF hexamer measuring the strain 
(er 5 ) dependence of the excitation gap Ai, Kr 5 the (q = 
0) inter- AF-hexamer interaction, and xr 5 (^) the strain 
susceptibility of a single AF hexamer. Here, we ignore 
the T-dependence in the background C44 (the hardening 
with decreasing T)22, because the nonmonotonic varia- 
tion in C±±(T) is very large compared with the harden- 
ing on cooling in C44. Fits of the experimental data for 
MgCr 2 4 and ZnCr 2 4 to Eq. §5$) are depicted in Fig. 3 
as solid and dotted curves, respectively. Here, the value 
of N = 3.45 x 10 27 m -3 is fixed in a first approximation^ 2 -. 
With the fitted parameters listed in the table beneath the 
figure, fits of Eq. (j5j) are in excellent agreement with the 
experimental data, reproducing the characteristic min- 
imum in C AA {T). For MgCr 2 4 (ZnCr 2 4 ), the fitted 
gap value of Ai = 39 K (34 K) in the PM phase is rather 
smaller than the gap of A ~ 50 K measured in the AF 
phase 3,11 , which implies that the suppression of spin fluc- 
tuations by the spin JT distortion makes the gap in the 
AF phase, A, rather larger than that in the PM phase, 
Ai. The negative value of Kr 5 = -19 K means that the 
inter- AF-hexamer interaction is antiferrodistortive. The 
coupling constant G 2 ,r 5 = \dA 2 /de T5 \ = 10200 K (9290 
K) is about three times larger than Gi 5 r 5 = 3600 K (3160 
K), indicating that the higher excitations A 2 = 136 K 
(111 K) couple to the trigonal lattice deformation more 
strongly than the lowest excitations Ai. 

According to Eq. (j2j), not only C44 but also C\\ should 
couple to the AF hexamers via the exchange striction 
mechanism, while Ct should be inactive. However, the 
nonmonotonic T dependence is invisible in Cu(T). In- 
stead, C\\(T) exhibit Curie- type softening due to the 
dynamical spin JT effect. We notice here that the mag- 
nitude of the nonomonotonic T dependence in C±±(T), 



especially the relative change in hardening on cooling be- 
low - 50 K down to T N in C 44 (T) (AC44/C44 ~ 1%), 
is rather smaller than the magnitude of Curie-type soft- 
ening in Cn(T) (ACn/Cii - 10%). Thus, in Cn(T), 
the weak nonmonotonic T dependence would be hidden 
behind the huge Curie-type softening. 



C. Hardening in AF phase 

As described above, the elastic anomalies in MgCr 2 04 
and ZnCr 2 04 in the PM phase strongly suggest the coex- 
istence of dynamical spin JT effect and dynamical molec- 
ular spin state, both of which strongly couple to the acous- 
tic phonons. Conversely, in the AF phase, C t (T) and 
C^T) exhibit the usual hardening on cooling without 
anomaly, as shown in Figs. 2 and 3, respectively 2 ^. Ac- 
cording to the INS experiments, while the spin JT distor- 
tion below T/v partially releases frustration, the dynami- 
cal molecular spin state persists even in the AF phase in- 
dicating the survival of some frustration 3 -^. Thus, com- 
bined with the above discussion on the elastic anomalies 
in the PM phase, we attribute the absence of anomalies 
in Ct(T) and C±±(T) in the AF phase to different origins: 
respectively, the spin JT distortion (generating the par- 
tial release of frustration) for Ct(T), and the decoupling 
of the acoustic phonons from the dynamical molecular 
spin state (producing the decoupling from the remaining 
frustration) for C±±(T). 



V. SUMMARY 

To summarize, ultrasound velocity measurements 
of MgCr 2 04 revealed elastic- mode-dependent sound- 
velocity anomalies in the PM phase. These elastic 
anomalies can be attributed to a coexistence of the dy- 
namical spin JT effect and the dynamical molecular spin 
state in the PM phase, the presence of which is compat- 
ible with the coexistence of the magnetostructural order 
and the dynamical molecular spin state in the AF phase. 
Further experimental and theoretical work is necessary 
to clarify the coexistence mechanism for the two different 
types of frustration effects in the geometrically-frustrated 
chromite spinels. 
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